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Abstract. In this note we discuss explicitly the structure of two simple set 
of zeros which are associated with the mean staircase emerging from the zeta 
function and we specify a solution using the Lambert W function. The argu- 
ment of it may then be set equal to a special N X N classical matrix (for every 
N) related to the Hamiltonian of the Mchta-Dyson model. In this way we 
specify a function of an hermitean operator whose eigenvalues are the "trivial 
zeros" on the critical line. The first set of trivial zeros is defined by the rela- 
tions Im(C (| H-i-t)) =0ARe(C(| +i-t)) ^ and viceversa for the second 
set. (To distinguish from the usual trivial zeros s = p + i ■ t = —2n, n 1 
integer) 

This (heuristic, non rigorous) research note is dedicated to the international 
Swiss -Italian mathematician and physicist Professor Dr. Sergio Albeverio on the 
occasion of his seventieth birthday; a friend and for years the scientific director 
of Cerfim (Research Center for Mathematics and Physics of Locarno), situated 
opposite the "RivellinoTA 



1. Introduction: a search for an hermitean operator associated 
with the rlemann zeta function 

There is much interest in understanding the complexity related to the Riemann 
Hypothesis and concerned with the location and the structure of the non trivial 
zeros of the Riemann zeta function £(s) where s = p + i ■ t is the complex variable. 
Following a suggestion of Hilbert and Polya, in recent years many efforts have been 
devoted to a possible construction of an hermitean operator having as eigenvalues 
the imaginary parts t n of the non trivial zeros of £ (£ being meromorphic, the zeros 
are countable). These are given by the solutions of the equation ( (p n + i ■ t n ) = 0, 
n = 1,2, . . . ,. If p n = \ for all n, then all the zeros lie on the critical line (the 
Riemann Hypothesis is true) ; the program is then to find an hermitean "operator" 
T such that T ■ tp n = t n ■ ip n in some appropriate space (ip n would be the n th 
eigenvector of T) . There are today many strategies in the direction of constructing 
such an operator and in the sequel we will comment on some (among many others) 
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very stimulating works on the subject. In pQ, Pitkanen's heuristic work goes in the 
direction of constructing orthogonality relations between eigenfunctions of a non 
hermitean operator related to the superconformal symmetries; a different operator 
than the one just mentioned has also been proposed in [5] by Castro, Granik and 
Mahecha in terms of the Jacobi Teta series and an orthogonal relation among its 
eigenfunctions has also been found. In the rigorous work by Elizalde et.al [3] some 
problems with those approches have been pointed out. In a work of some years ago 
Julia [1] proposed a fermionic version of the zeta function which should be related 
to the partition function of a system of p-adic oscillators in thermal equilibrium. In 
two others pioneering works of these years, Berry and Keating E] proposed an 
interesting heuristic operator to study the energy levels t n (the imaginary parts of 
the non trivial zeros of the zeta function). The proposed Hamiltonian has a very 
simple form given, on a dense domain, by: H = p ■ x + §, where 



(1) P = 



1\ d_ 
ch- 



in one dimension. As explained by the authors, the difficulty is then to define 
appropriate spaces and boundary conditions to properly determine p and H as her- 
mitean operators. In such an approach the heuristic appearance of "instantons" is 
also discussed. In another important work Bump et al. [7] introduced a local Ric- 
mann Hypothesis and proved in particular that the Mellin transform of the Hermite 
polynomials (associated with the usual quantum mechanical harmonic oscillator) 
contain as a factor a polynomial p n (s), corresponding to the n-energy eigenstate 
of the oscillator, whose zeros are exactly located on the critical line a = \. The 
relation of the polynomials p n {s) with some truncated approximation of the en- 
tire funcion £(s) (the Xi function), related to the Ricmann zeta function seems to 
be still lacking. Others important mathematical results concerning the non trivial 
Riemann zeros, have been obtained by many leading specialist (see among others 
the work by Connes [8] and the work by Albeverio and Cebulla j9]). For a recent 
work on the xp hamiltonian see G. Sierra [TU] , 

Let us also mention that for the nontrivial zeros of zeta an interesting equation 
has been proposed originally by Berry and Keating in [5]. In fact, remembering the 
definition of £(s) where £(s) is the Riemann zeta function, given by: 

(2) £(«) = i.«.(*-i).7r~ f .r(|).c(«) 

the equation for possible zeros of £ proposed in [5] is given by: 

S 1-5 

7T2 7T 2 

(3) —^x + —n — c = 

r(f) r(i^) 

As stated by the authors, Eq. ([3| could be considered as a "quantization condition" . 
Unfortunately, as mentioned in[5], Eq. ^ possesses complex zeros and so can not 
be used to provide an hermitean operator which would generate the non trivial 
zeros of (. The content of our note is concerned with the "mean staircase" of the 
Riemann zeros: we first specify the two set of trivial zeros on the critical line related 
to it and point out an explicit construction using the Lambert W function; then we 
introduce a specific argument (anxn hermitean matrix H , describing a discrete 
harmonic oscillator with creation and annihilation "operators" a and a* such that 
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[a, a*] = —2) into the Lambert W function. We obtain then, for the trivial zeros, 
the goal that the "Polya-Hilbert program" has for the non trivial zeros. 

2. The mean staircaise of the Riemann zeros and the trivial zeros on 
the critical line associated with it 

Let £(s) be the xi function given by If N(t) denotes the number of zeros of £ 
in the critical strip of height smaller or equal to t, and if S(t) = — arg K (| + it)), 
then [TT] 



(4) N(t) = (N(t)) + S(t)+0 
where 

(5) W)> - i • (* (£) " l) + I 

(N(t)), the "bulk contribution" to N, is called the "mean staircase of the zeros" 
(cfr. The fluctuations of the number of zeros around the mean staircase, are 

given by the function S(t). It is known [IT] that S(t) = O(lni) without assuming 

RH while, assuming RH is true, it is known that S(t) — O ( dj^j^ ■ At this point, 
since our remark has mainly to do with (N(t)), we will set S(t) = in Eq.Q. 
We shall study the relation N(t) = (N(t)) + 0{\). The two sets (which we call 
here "trivial zeros on the critical line" ) of interest are defined by the above mean 
staircase as follows. The first set is given by the zeros of Im£ (| + it) alone i.e. 
such that Re M | + it) ^ (the first set of trivial zeros on the critical line) . The 
second set is given by the zeros of Re£ (| + it), such that Im£ (~ + it) ^ 0. For 
the first set: 

(6) ImcQ+<) =0ARecQ+<) 7^0 

Then since Im£ (| + it* n ) = we have that 7r _1 arg^ (| + it%) = ir^ 1 (—nri) and 
7T _1 arg£ = is given by those t* n such that 

(7) N(t* n ) - (N(t*J) = § • (in - l) + I = n, ninteger > 1 

at large values of t or n. So the nonlinear equations to be solved which should give 
the values where only Im£ (| + it) vanishes, i.e {t* n }, is given by: 

(8) (N(t* n )) = n, ninteger ^ 1 
while for the second set 

(9) (N(t* n *)) = n- i ninteger ^ 1 

(The first set has been known for a long time and constitutes the Gram points, 
sin(6>) = 0, where 9 is the phase of the ( function, while for the second set one has 
cos(0) = 0). 

We note, the values of interest are given by the abscissa of the intersection 
points between the staircase (Eq.(j5|) and the two functions 7r _1 arg£ Q + it) an d 
7T _1 arg £ + it) — | . The plot of Fig 1 illustrates the situation for some low lying 
zeros. The values for £* lie mostly in between the exact value of the Riemann zeros 
i n _i and t n , but it is known that the Gram law fails for the first time at t = 282.4 
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( "first istanton" according to j5] ) . The solution of the above equation which gives 
t*, £** using a very special function (the Lambert W function, see [T^]) is given 
below. 



5- 




10 20 



-1- 



Figure 1. The plot of (N(t)) (continuos curve), of N(t) (full stair) 
and N(t) — \ (intermittent stair) 



3. An exact solution for the sequence i* and t* 
The equation corresponding to (|8|, may be written in the form 



(10) 



( — 

and the equation corresponding to ^ in the form 



(id 



\2Tie) 



so that introducing the new variables x — exp 8 J resp. x — exp ^— — | — M we 
obtain the equation (from (10 1 and (11), x > 0) 
(12) W(x) ■ exp(W(x)) = x 

The function W(x) is called the Lambert W function and has been studied exten- 
sively in these recent years. In fact such an equation appears in many fields of 
science. In particular the use of such a function has appeared in the study of the 
wave equation in the double-well Dirac delta function model or in the solution of 
a jet fuel problem. See [T^] for an important work on the subject. Moreover the 
Lambert W function appears also in combinatorics as the generating function of 
trees and as explained in [T2] the W function has many applications, even if the 
presence of the W function often goes unrecognized. 

The Lambert W function has many complex branches; of interest here is the 
principal branch of W which is analytic at x = 0. So, the solution of (pi |9|) is given 
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by 

(13) t* n = 2ne ■ exp (w 

(14) t~ = 2ire ■ exp (w 

We have thus specified, with the help of the LambertW function, the sequences {t* } 
resp. {t'n*}, which are the zeros of Im (C (\ + i ■ t)) such that Re (C (| + % ■ t)) ^ 
resp. Re (C (3 + i ■ t)) = and Im (C (| + i ■ t)) ^ 0. 

It should be noted here that in Eq.(J8|, n, which would correspond to the exact 
value of a true zero value t n (non trivial zero) of the £ function would not be an 
integer n or n — \ since we have replaced in Eq. ^ S(t) by zero. For the first few 
low zeros (the true zeros), it may be observed numerically that the corresponding 
values, let say n* , are randomly distribued mostly between two consecutive integers, 
but the mean values are nearby the integers plus \. A calculation with some zeros 
gives a mean value of 0.49 instead of 0.5. So, in average it seems that the behavior 
of the true zeroes t n "follows" more the pattern of the set t* n * . In the similar way 
the zeros of the first set, i.e. f*, lie mostly in between two non trivial zeros of C but 
of course it is known that there are very complicated phenomena associated with 
the chaotic behavior of the non trivial zeros of the Riemann £ function. 

As an example, the first of the istantons, corresponding to n = 126, cited above, 
is located at the value of t = 282.4 .... On the table below we give the values of i* 
and of t n of a true zero around t = 280. 

t 126 = 279.22 9 25 
t* 126 = 280.80 2 46 
t* 127 = 282.4547596 
*12 7 = 282.4651147 
t U8 = 283.211185 
t* 12S = 284.1045158 
ti29 = 284.8359639 
Table 1. 





From those numerical computations we see that two consecutive zeros of Im(£) 
alone arc followed by two consecutive true zeros, that is t\ 2 7 anticipates £127- 
The difference between the two subsequent t values is very small and given by 
At = 0.0103. The phase change is given by iir as illustrated on the plot of 
Im (in (£ (3 + it))) (step curve) and that of Im (£ (| + it)). 

For the first 500 energy levels, that is for values of t from to t = 811.184. . . 
(level number n = 500), it may be seen that there are 13 istantons (in the language 
of [5]), all with a Maslov phase change of +iir or of — iir. The width is usally small 
but it is larger for the istanton located at t = 650.66 (n corresponding to 379), 
where this time At — 0.31 .... Returning now to the trivial zeros (the two sets i*, 
<** defined above), we note the element ar relation which follow from |8]) and ([9]), 
and given by: 

t* + tt , -, 
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Figure 2. The first istanton 



and 



(16) 



r* i + r , ! 



= t* 



Eq.(Jl5 16 1 say that the zeros of the real part alone are obtained by those of the 
imaginary part alone by simple average and viceversa. The two sequences are 
regularly spaced and the mean distance between two trivial zeros at the height t, 
as the mean staircase indicates (Eq.([5])), is given approximatively by: 

(17) ' 2n ^ 



(N(t)) log i 



2tt ) 



log(rc) 



for t and n large 

Before proposing an hermitean operator for the sequences of the trivial zeros it 
is important to investigate a possible "quantization condition" for the non trivial 
zeros. For this we start with the Riemann symmetry of the £ function. 

From the exact relation for the ^-function given by: 

1 



(18) ^(s) 



1 



= ^1 - S) = -TP" 



s G CI, we have that 

(19) 71-- * r 



(l)cw 



T 



C(i-*)(i-*)(i-*-i) 



C(i-s) 



In equation (19) we limit ourselves to consider the values s = p+it=^{ +it, 
and thus 1 — s = \ \ —it; moreover we are interested in high values of t so 
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that we may use the Stirling's formula for the Gamma function given by: 
(20) T(x) = (2tt)5 x x -h~ x 



as x — > oo. From (19 1 and (20 1 we then obtain (asymptotically for t — > oo) 

it 



,2-1, exp(^((^ )Ml 



1 11 2!) ( ln (i 



I j - r ] +iarg ^ 
I I - I ) +zarg U 



it 



Since 



exp ( iarg [ C ( - T ~H 



exp [ia,rg[C[~i+it]+m 



exp ( i arg [ C ( - j 



we then have, taking the limit p = 5 |= 5, that: 



(22) 

Thus W 

hence 
(23) 



cos(*) = where* = ^ (in ( — - 
2 V V 2n 



axgfcf-+^ 



7r (n + |) . We then obtain: 



h ( ln 



27 ( ln \h 



1 1 



arg C 



- 1 



7 1 



arg C 



it \ = n 



it) = n 



Eq. ( |22[ ) may be seen as an approximate "quantum condition" for the true Ricmann 
zeros, but it is only a consequence of the Riemann symmetry (Eq. 19 1. In fact, if 



in equation ( |22[ ) we neglect the last term arg(£), then (22 1 has as a solution the 
second set of trivial zeros {t**}. It is true, as remarked by Berry and Keating, 
that their Eq.(|3| has complex zeros wich are not the Riemann zeros, but it should 
be remarked that if in Eq.|3| we set Re(s) = ~ then Eq.|3]) reduces to Eq.(p2) 
without the fluctuation term arg(C); so the solution of Berry and Keating Eq?l3) 
for Re(s) = | is the same as the second set of trivial zeros {t^*} we have specified. 



Below the plots of the left hand side of Eq.(22), with and without the term 
argC(| + it). As an illustration, we may observe on the plot the first istanton 
discussed above and the second one. In fact the maximum of the function which 
gives t* (Eq. [22] without the term arg () is outside the plot of the step function 
given by (22 1 (the true function) . This is visible on the plot near t = 282 and near 
t = 296 (the second istanton). This conclude our remark on ^ and Eq.(22l. In 
the next section, we shall costruct an hermitean operator whose eigenvalues are the 
trivial zeros of the zeta function on the critical line. 



Now in (13) and (14 1 the value of a trivial zero (t* or £**) is given through his 



index n by means of the Lambert W function so that such zeros are related in a 
non linear way to the integers n, i.e. in principle to the spectrum of an harmonic 
oscillator. So, for the trivial zeros, no boundary condition is needed here, since 



they are obtained by means of (13 1 and (14) in the large t limit. At this moment 



we are free to introduce a hermitean matrix which may generates the trivial zeros. 
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Figure 3. Plot of the function cos(^) of Eq.(22l with the term 
arg (C (| + it)) (step function) and without that term. 



4. An Hermitean operator (matrix) associated with the mean 
staircase (trivial zeros) of the rlemann zeta function 



As remarked above, in (13 1 and (14 1 the only "quantal number" is the index 



n of the trivial zeros and the construction may be given using a hermitean n x n 
matrix H, for any n, at our disposal and related to the classical one dimentional 
many body system whose fluctuation spectrum around the equilibrium positions is 
that of the harmonic oscillator. In fact, the one dimensional Mehta-Dyson model 
of random matrices (which may be seen as a classical Coulomb system with n 
particles) has, at low temperature an energy fluctuation spectrum given by the 
integers and it is possible to introduce classical annihilation and creation operators, 
as studied in |13j (a short discussion is presented in the Appendix). The matrix 
elements of the associated hermitean matrix are then functions of the zeros of the 
Hermite polynomials; in this case we do not have a Hilbert space and no Schrodinger 
Equation will be associated with the Lambert W function. Another direction, i.e. 
that of introducing a Schrodinger Equation to describe the trivial zeros may in 
principle be obtained as an application of the results given by G. Nash [2]; this 



because for large n, as it is known, (13 1 and (14 1 give the behavior pag. 214) 

related to the asymptotic behavior of the LambertW function: 



(24) 



tn 



27rn 
ln(n) ' 



and thus the spectrum appears in fact as a one where the associated Schrodinger 
Equation contains a Gaussian type of potential [2]. Here we will consider the 
matrix formulation: the point may seem to be somewhat artificial but the hermitean 
matrix we will use (specified in the Appendix) is related to the Mehta-Dyson model, 
the "starting point" of the random matrix theory. To do this, we begin to write 
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(10 1 in a slighly different form using the Stirling formula for the Gamma function 



of real argument given by: 

T(x) = (2ir) 2 x x ie ^asx-^oo 
We then have that, as t — > oo, 

f ' • l \) ' 'J' ' ^ 7 1 ^ 7 



(25) In T + - = In 1 + - + - In (2tt) - 

1 ; V \2ttc 2JJ 2ttc \2itc J 8 2 1 ' - 

1 7 

= n* + - In (27r) = n* + 6 

2 8 

where = \ In (2tt) - \ . 

Thus introducing the operator T — T(H) whose eigenvalues should be the trivial 



energy levels (for the first as well as for the second set defined by (13 1 and (14) as 
well as H, the hermitean matrix given in the Appendix and related to the Mehta- 
Dyson model, we may write the following heuristic matrix equation: 

(26) r(£ + 5)= e " + " 



where / is the unit matrix. Eq.(26l is the equation for T, giving the trivial zeros. 



The inversion of this formula (if it is possible to take it) yields hcuristically: 
(27) T = T(H) = 2tt (V 1 (e H+e ) - Q 

To conclude, if Hip n = (n + 1) <p n , where ip n is the n th eigenfunction of H, then 



(28) Tip, 



= 2tt (V- 1 (c H+e ) Vn = 2n (t- 1 (c"+ 9 ) - Vn = t nVn 



where t n = t* n resp. i**. (in thats latter case with 9 lowered by 1/2) 

Of course Eq.(26 1 for the operator T is more appealing than (13] 14| (where n is 



replaced by H and t* resp. i** are replaced by T) due to the combinatorial nature 
of the Gamma function, but the eigenvalues of the operators are the same in the 
"termodynamic limit", t — > oo. 

Remark: If one consider the usual map z —t 1 — ^ then the critical line s = \+it 
(t £ R) is mapped onto the unit circle \z\ = 1; the two sets of trivial zeros {i* } and 
{£**} have as accumulation point z — 1 (as n — > oo), which is the same accumulation 
point for the real zeros of the £ function given by z n — 1 — — 1 + ^, as n — > oo 
(see Fig Q. 

Neglecting the real zeros {z n }, Fig. |4] illustrate by means of two sets of trivial 
zeros {£*} and {t^*} the Lee- Yang Theorem for the zeros of the partition function 
for some general spin lattice system studied in statistical mechanics. If RH is true, 
then all non trivial zeros of £(s) shall be located at the same circle \z\ — 1, with 
z = 1 as accumulation point. 

5. Conclusion 

In this note we have specified an operator equation for the operator T, having as 



eigenvalues the trivial zeros of the £ function on the critical line (13 and 14 1; such 
zeros are the so called trivial zeros (mean staircases), and are related in a strong 
non linear way to the eigenvalues of a discrete harmonic oscillator described by H 
in the large t limit. The introduction of H may seem to be artificial: nevertheless 
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kl = i 




Figure 4. z- plane 



the two sets of trivial zeros ({£*} and {t^*}) are the eigenvalues of T(H). In a 
subsequent note we will present a study of another sequence of zeros possibly more 
connected with the true Riemann zeros. 

Update 

Very recently G. Sierra and P.K. Townsend |15j introduced and studied an in- 
teresting physical model (a charged particle in the plane in presence of an electrical 
and a magnetic potential). In particular, the lowest Landau level is connected with 
the smoothed counting function that gives the average number of zeros, i.e. the 
staircase which here we have studied, by means of a classical one-dimcntional model 
of N interacting charged particles. 

Note added 

Very recently Schumayer et al [TB] constructed (in particular) the Quantum 
mechanical potential for £(s) zeros, with the first 200 energy eigenvalues (non trivial 
zeros) . It is expected that the same form of a quantum mechanical potential would 
appear using only the two sets of zeroes we have discussed in this note. For the 
construction of an Hamiltonian whose spectrum coincides with the primes, see also 
the recent work of S. Sekatskii [TT] . 
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6. Appendix 

The hermitean matrix associated with the Mehta-Dyson model, H: discrete 
annihilation and creation operators associated to H whose spectrum is given by 
the set of integers (1,2, ... ,n), for any finite n. 

In Ref [13j [18] it was studied the one dimentional Mehta-Dyson model defined by 

N ( \ 

the Hamiltonian E = \ \yl — Yl 1°8 {\Vi ~Vj\) ) > where yi is the position of 

i=0 \ i<j<N J 

the i-th particle on the line; then the fluctuation around the equilibrium positions 
(these are given by the zeros of the Hermite polynomials of degree N, where N is 
the number of particles on the line, for every finite N), i.e. the harmonic fluctuation 
spectrum is given by the eigenvalues of the hermitean N x N real matrix whose 
elements are given by: 

Hi} = l^-^-l 2 * ^ 3 



< 



Hvi = l + J^^F 1 = j 



i,j = 1,.. ., N, where now the Xi arc the "equilibrium positions" i.e. the zeros of 
the Hermite polynomials of degree N. 

The spectrum of H is given exactly by the integers (1, 2, ... N) for every finite 
N and the eigenfunctions are given in terms of the Mehta-Dyson polynomials of 
order 1 up to N . The Hamiltonian describing the harmonic fluctuations takes then 
the form [13] : 

H = N ■ I - -aa* 
2 

where / is the unit matrix of order N and a, resp a*, are the discrete annihilation 
and creation operators (matrices of order N x N) which satisfy the commutator 
relation [a, a*] = —2. 

Moreover [H, a*] — a* and [H,a] = —a. 

If Xk is the fc-ten eigenvector of H with eigenvalue the integer k, one has: 

a Xk+\ = Xk+2 

and 

aX k+1 = 2{N-k) X k 

Explicitly, if X k = (ipik [x\) v?Nk ( x n)) is the fc-ten eigenvector, where <fk(x) 
is the k-th Mehta-Dyson polynomial of argument x, then 

i s d ^<Pk+i (xi) - ifik+i (xi) 
a<Pk+i (xi) = -t— [tpk+i (xi)) = 2^ 7 — n 

CLXi . I X± Xi I 

and 

a*(p k+ i (xi) = ^2xi - <fk+i (xi) 

a and a* are the two discrete annihilation and creation operators of H. 

H as above with TV = n may be used to give the first n trivial zeros of the first 
set in (13 1 i.e. t\ . . . i* while H + \ may be used for obtaining the first n trivial 



zeros ofthe second set in ( 14 1 i.e. t** . . . t** in the discussion on the mean staircases 
given in Section 4 above. 
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